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ABSTRACT 

The tuning job of the compensator for levitating flexible rotors supported by active magnetic 
bearings (AMB) concerns providing a good damping effect to the critical speed modes while avoid- 
ing the spillover problem on the instability of higher bending modes. In this paper, an idea for 
design of the control law of the compensator based on utilizing a third order low pass filter (LPF) is 
proposed to essentially enable elimination of the spillover instability. According to the proposed 
design method, good damping effects for the critical speeds are obtained by the usual phase lead/lag 
function. Stabilization for all of higher bending modes is completed by the additional function of 
the 3 rd order LPF due to its phase lag approaching about -270 degrees in the high frequency domain. 
I his idea is made clear by experiments and simulations. 


INTRODUCTION 

To design control network driving active magnetic bearings (AMB) for flexible rotor levitation it 
is important to consider rigid mode levitation, controllability of flexible modes at critical speeds 
and a stability margin for high frequency bending modes. For instance, a design of super critical 
compressors supported by AMBs illustrates this point. In the first design phase, rigid mode levita- 
tion performance is determined. Next, the well damping ratio of the first bending mode necessary to 
pass the critical speed is determined according to O-value criteria. Finally, the PID controller 
transfer function is tuned so as to satisfy these two requirements and so as to avoid higher frequency 
instability beyond 2nd bending modes. It is called a spillover problem. 111 

Usually, the PID control law with optional notch filters and/or low pass filters is used for tuning. 
Otherwise, the increase of the internal damping of the rotor is improved. However, this design 
method based upon the phase lead/lag function is not free from the spillover problem. However the 
presented method can provide a highly robust control system for the higher modes and satisfactory 
performance for the controlled modes. It is possible to skillfully combine the commonly used 

phase lead/lag function with the 3rd order LPF. How to design the 3rd order LPF is a main point of 
this paper. K 

The phase lead/lag function can possibly provide enough damping performance to the rigid and 
first bending modes. It can be done by placing these critical speeds in the phase lead domain. If a 
flexible rotor system requires higher damping than one provided by the phase lead function, some 
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optional function can be used just for passing the bending critical speeds According to our expe- 
riences, an optional network, called N-CROSS, is recommended.'-' It works only around the critical 
speeds of the shaft bending modes during the rotor rotation. Then, the problem is how to design 
the compensator featuring no spillover to potentially stabilize all the eigen modes. 

The proposed design method features the addition of the 2nd order LPF to the usual PID compe 
sator in order to realize this 3rd order LPF. The center frequency of the 2nd order LPF coincides 
with the eigen frequency of the rotor restricted by pin-pin boundary at AMB portions. This eigen 
frequency is equal to the anti-resonance frequency located between the ist and ~nd bendmg^modes 
Since the total function of the controller network becomes a 3rd order LPF the phase lag forwards 
to -270 degrees, i.e., a 90 degree phase lead, in the high frequency domain. Therefore, all the higher 
natural frequencies beyond the 2nd bending mode are essentially stable, i.e., completely free from 

This idea ^theoretically explained for a rotor system with one AMB support. The effectiveness 
of this design method is proven by numerical simulations and experimental results. The rotation 
test is done with a low enough Q-value to pass the rigid and first bending critical speed without he 
high frequency instability. The generalization of this design concept is finally proposed for the 
flexible rotor borne by two AMBs. 


A FLEXIBLE ROTOR WITH ONE ACTIVE MAGNETIC BEARING 


A flexible rotor shown in Fig.l is selected for explaining the fundamental idea of the design 
method of AMB compensator. The rotor is supported by a ball bearing at the right end and by an 
AMB at the other. The ball bearing placed here is represented by the pin boundary condition to 
explain the idea clearly. The problem then focuses on to how to design the compensator driving the 

AMB at the left side. 



rotor weight = 18 DN 
Free-free natural frequency = 

47 Hz (1st) 175 Hz (2nd) 

Rated speed = 5,000 rpm 


A gap sensor placed at the left side detects the rotor radial displacement which feeds its signal to 
the controller. The signal from the controller drives the magnetic bearing force through the power 
amplifier having a transfer function of unity. The rotor is connected to a motor with a flexible 
coupling. The flexible coupling is flexible enough to satisfy the free boundary condition at the right 
end of the rotor as shown Fig.2(a). Gyroscopic effects are neglected here. 

The displacement vector is noted by X=[X \X,j: X 0 for the rotor displacement at the AMB portion 
and X for the displacements of all portions of the rotor except the AMB portion. The former is 
called 'the boundary displacement and the latter called the inner displacement vector. In the AMB 
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servo-eontrol system, the boundary displacement X is measured and the AMB force U acts upon 
the rotor according to a control law. Then, the equation of motion of the rotor-bearing system can 
be written in matrix form as follows: 
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where M. and K.. (i,j=l,2) indicate mass and stiffness matrices of the rotor system, respectively 
noted by 1 for the inner portion and by 2 for the boundary AMB portion. 

Any damping factors, e.g., shaft material damping, mechanical dampers, are neglected In the 
authors opinion, the consideration of such system damping factors, in addition, will potentially 
estroy the realization of practical robust control. Rotordynamics designers want to consider the 
compensator design for the AMB, similarly, the design for oil-film bearings, i.e., neglecting anv 
system damping. b 6 3 


(c) Reduced Model 



Deflection Mode 




(b) Transformation modes 
Fig. 2 Modelling by quasi-modal transformation 


SYSTEM MODELLING 

The equation of motion (1) of the rotor is formulated in FEM matrix forms on meshing by the 
shaft beam element. The quasi-modal system reduction in a category of the modal synthesis is 
applied to the original system.'- 1 The quasi-modal transformation modes are determined by two 
kinds of modes as shown in Fig.2(b): y 

(1) Pure bending mode 0 obtained by eigen mode of the rotor if the pin condition is imposed 

. . at the AMB portion. 

(2) Rigid mode 6 obtained by deflection mode if the rotor is lifted to give the unit displacement 

at the AMB portion. ' 
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The quasi-modal transformation is defined by the following equation: 
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where the variable s is a weighting value indicating the magnitude of the pure bending mode ob- 
tained in the rotor bending critical speed mode. 

From the result of the transformation, the following equation of motion of the reduced model 
having 2 degrees of freedom is completed: 


m* 


s 

+ 

■ k* 


S 


" o ' 







i . 

m c - 


. *2 J 


0 . 


- X 2 _ 



where ■* = *%*. n«=*>Mi<5, m^=M 2 + <5 l M : c5 , k^cop 2 and 

co? = the natural frequency of the pure bending mode under the pin 
boundary condition at the AMB portion. 


This equation of motion having the mass interaction between the pure bending mode and the rigid 
mode is equivalently converted to the spring interaction. This conversion is possible to provide the 
reduced vibrating model visually composing the mass-spring connection, if a new absolute dis- 
placement x* is introduced as s=a(x* - x 2 ): 
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The obtained model is illustrated in Fig.2(c). The masses of meq 6 and meq indicate the equiva- 
lent masses of the rigid mode and the pure bending mode, respectively. 


CONVENTIONAL DESIGN CONCEPT FOR COMPENSATOR 

The negative spring effect exhibited by the AMB is neglected. It can be easily compensated by 
certain gain of the proportional action. The equation of motion (5) of the model is rewritten by the 
state equation form as follows: 

X = AX + bu 
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The controller design based upon 
the modem control theory combin- 
ing the full order observer and the 
optimum feedback gain due to 
LQR is applied to the model in the 
manner as shown in Fig. 3. [4 l 
One example of the transfer func- 
tion of the compensator itself is ob- 
tained as shown in Fig.4. The over- 
all configuration of the transfer 
function is fundamentally similar to 
PID action, but down slope gain of 
1st order LPF at the high frequency 
domain. 



The free-free eigen fre- 
quency is indicated by the 
bottom points of the dotted 
curve. Therefore, the damp- 
ing ratio of the rigid and 
first bending mode of the 
global system can be well 
provided. However, higher 
bending modes than the 2nd 
will be easily unstable due 
to the phase lag at this 
high frequency domain. 

In fact, the phase curve on 
the Bode diagram changes 
from the phase lead region 
to the phase lag region be- 
yond the peak portion of the 
gain curve. This phase lag 
forwarding to -90 degrees 
provides a negative damp- 
ing effect for high frequency 
eigen modes of the global 
system. This is the reason 
why the spillover instability 
is induced. 


Fig. 3 Conventional control layout 



Fig.4 Controller transfer funnction (Conventional design) 
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PHASE degree 





PROPOSED DESIGN CONCEPT FOR COMPENSATOR 

The idea of the new design comes from the fact that if this phase lag forwards to -270 degrees 
insTead of -90 degrees, the compensator can provide the positive damping to the system i even n the 
hieh frequency domain. 151 This idea can be realized by replacing the behavior of the 1st order LPF 
oHhe conventional compensator by one of the 3rd order LPF In other words, it is easily completed 
hv the cascade combination of the present type plus a 2nd order LPF as shown in Fig.5. 
b The type of the 2nd order LPF that must be added is now discussed. It is recommended I tha the 
center frequency of the 2nd order LPF is tuned with the anti-resonance frequency, i.e., the eigen 
frequencyof the^ure bending mode. This is because the rotor does not react to the resonance of the 

this 2nd order LPF, the required compensator is automatically designed LOnc of ^ teansfer f 
tions of the proposed compensator is shown in Fig.6 including the 2nd order LPF. I he dotted 
indicates the inverse function of rotor response against force input, called the rotor gain. The peak 
!^f the dotted curvehdicates the anti-resonance frequency in which the amplitudes of the corre- 

sDondine AMB portions are zero, i.e., nodes on the pin-pin eigen modes. , 

In this figure, the phase lead covers the natural frequency of the rigid and first be nding modes to 
provide the well damping effects at the critical speeds. The phase lag starts beyond 80 Hz, but it 
quickly passes -180 degrees before the 2nd bending mode and it finally moves forward to -27 
degrees. Therefore, positive damping is provided to all bending modes of the global sys em ey n 
the 2 nd mode, ensuring stability. As 
obviously shown in this phase curve, i o 0 

all of the higher frequency modes are 
definitely stable. There is no 
spillover problem. This controller 
will be robust from the practical 
viewpoint of the AMB compensator 
design. 

The result of the eigenvalue of the 
global system is shown in Table 1. 

Clearly all modes are completely sta- 
bilized even with no material damp- 
ing. 
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Fig. 5 Control layout 

(Proposed design) 


Fig.6 Controller transfer function 
(Proposed design) 








Table 1. Eigenvalue of the Flexible Rotor 


Mode 

Orpm eigen value (£%, Hz) 

Open loop 

Closed loop 

Rigid mode 
1st bending mode 
2nd bending mode 
3rd bending mode 
4th bending mode 

0 . 0 
0 , 47 

0, 175 
0 , 285 
0 , 374 

61.65038, 20 

16.24556, 48 

2.82522, 169 
0.05888, 285 
0.00009, 374 


EXPERIMENTAL RESULTS 

A control law is made by digital means with the following specifications: 

- Bilinear S-Z transformation used for digitalization, 

- 8 KHz for sampling frequency and 

- DSP board (type ADSP320 made by Chubu Denki) 

In the experiment, the integral action is added for statically levitating the rotor. A proportional 
action is also included in addition, for cancelling the negative spring effect in the AMB. 

For experimentally checking the stability of the developed control law of the 3rd order LPF, it is 
compared with the commonly used PID control law as shown in Fig.7. In the case of the PID 
control, depending on the increase of the feedback loop gain, the instability of the 5th bending mode 
appears when the natural frequency of the rigid mode reaches 13 Hz. However, in the case of the 
3rd LPF control, even if the rigid mode natural frequency is increased to 24 Hz, the bending mode 
instability does not appear. This proposed controller is thus very robust. 
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As shown in Fig. 8, the vibration 
response curves are measured. 

Compared with the PID, enough 
damping effect is provided by the 
proposed controller for the rigid and 
first bending modes at the critical 
speeds. Small peak amplitudes are 
obtained compared with the con- 
ventional type. 

As suggested in Fig.9 of the rela- 
tion between the natural frequency 
versus the AMB stiffness, the com- 
pensator design concept is as fol- 
lows: 

- 1st eigen mode (rigid mode) and 2nd eigen mode (1st bending mode) can take enough damping 
ratio to pass the critical speeds owing to the phase lead. Consequently small Q-value design 
is completed. If necessary, the N-cross can be optionally added for making the Q-value of 
the 1st bending mode smaller. 

- The higher eigen modes than the 3rd (2nd bending mode) are stabilized due to the phase lag 
of about -270 degrees. The spillover problem completely disappears for any higher frequency 
mode. 

- The center frequency of the additional 2nd order LPF must be set with the anti-resonance 
frequency located between the 2nd and 3rd eigen modes (1 and 2nd bending modes). Even 
if the Q-value of the 2nd order LPF is high, the rotor does not react to such sensitivity as 
high as in the compensator network. 
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GENERALIZATION OF PROPOSED DESIGN CONCEPT 

The rotor borne by two AMBs at the left and right sides is quite common. The concept of the 
compensator design developed for one AMB type can be enhanced to two AMB types. The natural 
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BEARING STIFFNESS DN/na 2nd LPF 

Fig. 10 Natural frequency map (Two AMBs) 

frequency map is shown in Fig. 10. On 
the left side of the map the natural fre- 
quency of the rotor subjected to free- 
free boundary condition is plotted. On 
the right side the natural frequency is 
that of the pin-pin boundary condition. 

The right side frequencies indicate anti- 
resonance frequencies. 

For the type of two AMB there are 
two types of controller layouts: central- 
ized and decentralized. The centralized 
layout employing the translation and 
tilting mode controls are suitable for the 
control law of the 3rd order LPF. The 
separation into both the modes and two 
independent compensators are prepared 
asshowninFig.il. 

According to this separation, two Fig. 1 1 Centralized controller layout 

types of the compensator are indepen- 
dently designed, based on the design concept developed for the one AMB control : 

- type 1 : translation mode compensator to control mainly the parallel rigid mode and the 1st 

bending mode of the rotor. 

- type 2 : tilting mode compensator to control mainly the conical rigid mode. (The 2nd bending 

mode of the rotor must be included in the control frequency band, if necessary.) 
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The obtained transfer functions of each compensator are shown in Fig. 12 and 13. The center fre- 
quency of the 2nd order LPF coincides with the 1st anti-resonance frequency of the rotor under the 
pin-pin boundary condition, as shown in the natural frequency map of Fig. 10. 

P The stability analysis of the global system is shown in Table 2. All of the bending modes are 
stabilized without spillover problems, even if no material damping. Clearly the completed regu 
tors are very robust. 



Fig. 12 Transfer function for 
translation mode controller 
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Fig. 13 Transfer function for 
tilting mode controller 


Table 2 Eigenvalues of the Rotor (2 AMBs) 


Mode 

Orpm eigen value (£%, Hz) 

Open loop 

Closed loop 

Translation mode 
Tilting mode 
1st bending mode 
2nd bending mode 
3rd bending mode 
4th bending mode 

0 , 0 
0 . 0 
0, 154 
0, 331 
0, 555 
0, 786 

26.66492, 29 

18.48787, 31 

2.93966, 152 
0.00290. 331 
0.00091, 555 
0.00023, 786 
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CONCLUSION 


A design concept for AMB control based on a 3rd order LPF is presented. It is a key point that the 
center frequency of the 2nd order LPF within it must be set by the anti-resonance frequency of the 
rotor. The rotor system can be stabilized by positive damping due to phase lead for the rigid and 1st 
bending modes, and due to phase lag with about -270 degree for the higher bending modes. The 
completed compensator can potentially stabilize all of the bending modes with no material or struc- 

ure of damping. The effectiveness of this design concept is proven in the experiment using the 
rotor having one AMB. 6 

This concept is generalized for the flexible rotor borne by two AMBs, as this is a common type 
and its effectiveness is made clear by simulation. 
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APPENDIX 

, ? y u j ing T ^ ntro1 °hjector composing of the rotor and a 2nd order LPF, the compensator of the 
3rd order LPF is automatically designed. The corresponding state equation of this enlarged objec- 
tor as shown in Fig. 5 is written by the following formulas: 
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From the state eqation, a control law obtained by the LQR manner is thus shown in Fig.6. 
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